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Abstract 

We describe a technique that enables strong, coherent coupling between individual optical emit- 
ters and guided plasmon excitations in conducting nano-structures at optical frequencies. We show 
that under realistic conditions, optical emission can be almost entirely directed into the plasmon 
modes. As an example, we describe an application of this technique involving efficient generation of 
single photons on demand, in which the plasmon is efficiently out-coupled to a dielectric waveguide. 

PACS numbers: 32.80.-t, 03.67.-a, 42.50.Pq, 73.20.Mf 
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The fields of quantum computation and quantum information science have spurred sub- 
stantial interest in generating strong, coherent interactions between individual quantum 
emitters and photons. Such a mechanism would enable quantum information to be passed 
over long distances, which is essential for quantum communication and would facilitate 
the scalability of quantum computers |3| . The required coupling between emitters and pho- 
tons is difficult but has been achieved in a number of systems that reach the so-called "strong- 
coupling" re gil ue of cavity qU a„ tUm eiectrooyuarnics QQ. Receutiy, sever. approaches to 
reach this regime on a chip at microwave frequencies have been suggested (JO, 0, |l| and 
esperiruentaHy observed j^. A key feature of these approaches is the use of conductors to 
reduce the effective mode volume V e g for the photons, thereby achieving a substantial in- 
crease in the coupling strength gocl/VKff- Realization of analogous techniques with optical 
photons would open the door to many potential applications in quantum communication 
and in addition lead to smaller mode volumes and hence faster interaction times. 

This Letter describes a method that enables strong, coherent coupling between individ- 
ual emitters and electromagnetic excitations in conducting nano-structures at optical fre- 
quencies, via excitation of guided plasmons localized to nanoscale dimensions. The strong 
coupling is possible due to the small mode volume associated with this sub-wavelength con- 
finement. We first examine the simple case of a conducting nanowire, where the relevant 
physical mechanisms can be understood analytically, and then consider an optimized ge- 
ometry of a metallic nanotip. Because of dissipative losses in metals the plasmon modes 
themselves are not suitable as carriers of information over long distances. We show, however, 
that the plasmon excitation can be efficiently out-coupled through evanescent coupling with 
a nearby dielectric waveguide, as illustrated schematically in Fig. [T] This can be used, e.g., 
to create an efficient single photon source, or as part of an architecture to perform controlled 
interactions between distant qubits. We find that single-photon efficiencies exceeding 95% 
are possible using relatively simple implementations. 

Surface plasmons jl^ are bound, non-radiative electromagnetic excitations associated 
with charge density waves propagating along the surface of a conducting object. For a 
smooth, cylindrical nanowire, it is convenient to write the fields in terms of cylindrical coor- 
dinates, Ej(r) = E l j jm ( y kj±p)e im ^e lk ^ z , where j = 1,2 denotes the regions outside and inside 
the metal, respectively. The functions Ej m are determined by Maxwell's Equations along 
with appropriate boundary conditions and are given in [ijj. The modes are characterized by 
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the longitudinal component of the wave vector ku , which is related to the free-space wavevec- 
tor ko = Uq/c, dielectric permittivity 6j, and imaginary transverse wave vector kj± = inj± 
by €jky = k 2 — k 2 ± . is related to the inverse decay length of the field away from the 
surface. The dependence of k\\ for a conducting wire on its radius R is shown in Fig. |2Ja), 
for a few lowest order modes. Throughout this paper, all figures and numbers presented use 
the optical properties of silver {e% = —50 + 0.6z) [17[ at a vacuum wavelength Ao = 1 /im, 
and assume a surrounding dielectric e\ = 2. 

For a conducting nanowire (| -^/eHkoR-^l, Ree2 < 0), there exists one, fundamental TM 



mode 



12J, Q| with axial symmetry (m = 0), while all higher-order modes are cut off. For 



this mode, kn, Kj±&C-i/R, indicating that the phase velocity v p h(xuJoR is greatly reduced 
while the transverse mode area A e ^ oc R 2 is localized to a region on the order of the wire 
size. The proportionality constant C_i depends only on 6^2 and is given by 

£2_ 2 

ei ~ (7 -log 2 + log C^idY 1 1 ' 



where 7~0.577 is Euler's constant. This sub- wavelength guiding 0: 



14. 



)lasmons in metal 



nanowires has recently been observed in a number of experiments 

At optical frequencies a conductor has finite losses characterized by Im €2, resulting in 
dissipation of the plasmon given by Im hi. Due to the tighter localization, 

Re fcii/Im fcn decreases and approaches some non-zero value («140) as R— >0, as shown in 
the inset of Fig. |2fa) for the fundamental mode. Physically Re fey/Im k\\ is proportional to 
the number of plasmon wavelengths the plasmon will travel before decaying. 

We now consider the emission properties near a nanowire of an oscillating dipole, which 
physically can be formed by a single atom, a defect in a solid-state system, or any other sys- 
tem with a dipole- allowed transition. This dipole can generally lose its excitation radiatively 
by emitting a photon, non-radiatively through dissipation of currents induced by the dipole 
in the metal, or into the guided plasmon modes. It is well-known that the corresponding 



spontaneous emission rates can be obtained via classical calculations of the fields |18|. For 
sub-wavelength systems this calculation further simplifies because it is sufficient to consider 
the quasi-static (H = 0) field solutions Q], the derivation of which we outline here. Given 
a point charge source at r' outside the metal, we write the static potential outside in terms 
of the free-space potential and a reflected component, $ ( r ) r ') + ^r( r > r ')j while $2( r ; r ') 
gives the potential inside the wire. Note that the potential due to a dipole po at r' is easily 
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found from the point source potential via $di P (r, r') = (po • V) $(r, r'). Since 3> r ,2(r, r') are 
solutions to the Laplace equation, V 2 $ rj 2 = 0, we can expand $ r in an appropriate basis, 

<Ey(r, r') = - y^(2-5 mfi ) cos m((f) - (f)') / dh a m (h) cosh(z - z')K m (hp')I m (hp), (2) 

71 m=0 ^° 

where a m (h) are arbitrary coefficients, and K,I m are modified Bessel functions. A similar 
expansion holds for $2, with the replacements a m — *> (3 m and K m —*I m . Expanding $0 = 
(47reoei|r — r'|) _1 in a similar basis, algebraic equations for a m , f3 m result by requiring that 
$0 + = $2 at the boundary and that eE^ is continuous here. The solution for a m is 
found to be 

(h) = 1 (eg - eJI'Jh^UhR) 

m{ ' e x exUh^K'JhR) - e 2 K m {hR)I' m {hR) ' 1 1 

The radiative decay rate is determined by finding the dipole contribution ~ J^'* r 3 to 
$^jp , r = (p • V) $ r at large distances r, which for a nanowire is due to the m = 1 term 
in Eq. (J2J). Physically 5p corresponds to an induced dipole moment in the nanowire, and 
leads to a modified radiative decay rate r ra d cx |po + <5p| 2 . The non-radiative and plasmon 
decay rates can be calculated via r io4a ioclm(p ■ E), where E = — V$di P (r, r')| r=r / is the 
total field at the dipole location. In the limit that the distance d between emitter and 
wire edge approaches zero, one finds that E diverges due to a substantial contribution from 
the sum over m in Eq. (j2J). This term is proportional to Im €2 and can thus be identified 
with non-radiative decay. At the same time, Eq. (jSJ) exhibits a pole in a that corresponds 
to excitation of the fundamental plasmon mode, and whose contribution to E yields the 
plasmon decay rate. For a nanowire, the radiative and non-radiative decay rates for a dipole 
oriented along p are given by [3| 

e 2 - ei R 2 



Trad /To 
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e 2 + ei (R + d) 



(4) 



r n on-rad/r ~ — „ /2 Im ( — ^ ) , (5) 



where To is the decay rate in uniform dielectric e\ |2fll |. while the plasmon decay rate is 

Kl{K lL (R + d)) 
Fpl/ro = ° pl • (6) 

a p i is a complicated expression but depends only on ei 2 . 

The scalings of the various decay rates can be intuitively understood. Away from the 
plasmon resonance (ei + €2~0), r ra d varies slightly from Tq due to a small change in the ra- 
diative density of states near the nanowire, while the 1/d 3 dependence in r non _ rad reflects the 



dissipation of a divergent current induced in the nanowire by the near-field of the dipole. The 
1/R? scaling for T p i can be understood from Fermi's Golden Rule, T p i = 2ng 2 (r, u )D(u ), 
where g(r,u ) oc 1 / y/ A e s oc 1/R is the position-dependent coupling strength between emit- 
ter and plasmon modes and D{ojq) oc {dui / dk^ 1 oc {ojqR)^ 1 is the plasmon density of states 
on the wire. 

The position dependence of the decay rates results in an optimal emitter distance d Q 
for which the probability of decay into plasmons r p i/(r p i + V) is maximized, where V = 
Trad + Tnon-rad denotes the total "non-plasmon" decay rate. For typical parameters, d Q is on 
the order of several R as R— >0, such that the emitter sits within the localized plasmon field 
but is not too close to the wire that dissipation become dominant. In Fig. Efb) we plot the 
probability of emission into non-plasmon channels, 1 — T v \/ (V + T function of R 

when the optimal d is chosen. This optimized "error" rate decreases monotonically as R—>0 
and approaches a small number oc Im e/(Re e) 2 indicating that the efficiency of emission 
into plasmons is ultimately limited by dissipative losses. As R— >0, one can achieve effective 
Purcell factors of r p i/T'^5.2x 10 2 for a silver nanowire. 

The nanowire is a simple system that illustrates the relevant properties of dipole emission 
and plasmon propagation. One immediately sees, however, that the increase in coupling 
achieved by letting R— >0 is accompanied by a decrease in the plasmon propagation length, 
which limits coherent processes of interest. Such limits can be circumvented with simple 
design improvements, which we illustrate specifically for the case of a metallic nanotip, 
assumed to have a paraboloidal profile given by p(z) = ^pvz (z > 0). One expects a similar 
enhancement of plasmon emission due to the nanotip, yet the tip can quickly expand to 
larger radii where losses can be significantly reduced. As in the nanowire case, one can 
calculate the emission rates based on the quasi-static field solution of a dipole near a tip, 
which is exactly solvable by working in parabolic coordinates. For an emitter located on-axis 
at position z = — \d\ < and oriented along z, we find that 



Here, a pl is a constant that again only depends on ei j2 , and C_x is the solution to Eq. 






(9) 
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From these decay rates one finds an optimized Purcell factor T p \/r'^2.5 x 10 3 as v — >0, as 
shown in Fig. 12(b), for the same material parameters as the nanowire calculation. In the 
case of a nanotip, however, we are primarily interested in the probability that the plasmon 
propagates up to some final radius R. We estimate this quantity by making an eikonal 
approximation based on the nanowire solution [22]. In particular, we assume that the 
plasmons are completely emitted into the end of the tip (z = 0) at a rate T p i, while the rate 
at which the plasmon emission successfully propagates to some larger radius R(z) is given 
by 

/ MR) \ 
T pl (R) = r pl exp 1-2 J \mk\\{p{z))dz\ . (10) 



Here k\\(p) is the nanowire solution at radius p. In Fig.|2(b) we plot Pe(R) = 1 — f p i(i2)/(r'+ 
T p i) as a function of R, optimized over the emitter position and v. This quantity corresponds 
to an effective error probability in which the plasmon mode is either not excited or fails to 
successfully propagate to some final radius R. For koR > 0.05 the tip leads to significant 
improvement in efficiency compared to nanowires of the same R. To check the validity of 
these approximations we have performed numerical (electrodynamic) simulations of dipole 
emission near a nanotip using boundary element method [23], with the resulting numerically 
optimized Pe{R) plotted in Fig. E(b). It can be seen that the theory agrees well with the 
numerics. Some typical simulation results are shown in Fig. 0(a). Here we plot for different 
emitter positions the quantity |Re (E x H*)|/T tota i, which is proportional to the energy flux 
of the system normalized by the total power output of the emitter. It can clearly be seen 
that choosing the optimal position results in efficient excitation of the plasmons at final 
radius R, while other positions can result in primarily non-radiative or radiative decay. 

Because of losses, the plasmon modes are not suitable as carriers of information over long 
distances. However, one can evanescently couple the plasmons to dielectric waveguide modes, 
which can form an architecture for a device to generate single photons on demand. Noting 



that the concepts behind single p 



presented in detail elsewhere 



24. 



roton gene ration with a single emitter in a cavity have been 



25l l26l |27| , here we illustrate a potential novel realization of 
a single photon device, shown in Fig.QJa). In this architecture, a single, optically addressable 
emitter such as a quantum dot sits atop a nanowire, which co-propagates over some length 
L ex with a nearby dielectric waveguide. In order to maximize the transfer efficiency into 
the waveguide, the longitudinal wavevectors k\\ of the plasmon and waveguide should be 
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approximately matched, and L ex should be optimized. In practice, for a given wire radius 
R, the matching condition results in some optimization of parameters such as the waveguide 
size. A similar setup using a nanotip instead of a nanowire is illustrated in Fig. ^b). Here 
the nanotip maintains a paraboloidal profile up to some final radius R, at which point it 
becomes a straight nanowire supporting plasmon modes with well-defined kn that can be 
easily be mode-matched with the waveguide. 

The out-coupling and single photon efficiencies can be calculated using standard mode- 
coupled equations based on Lorentz reciprocity j^. For simplicity, we take the dielectric 
waveguide to be a cylindrical optical fiber, whose modes can be calculated analytically , 
and set the surrounding and fiber core permittivities to be e\ = 2, e c = 13. In Fig. 0(b), we 
plot the optimized efficiency P for single photon generation as a function of R, for both the 
nanowire and nanotip, based on the decay probabilities obtained above and coupled-mode 
theory. We also include the predicted efficiencies using the boundary element method sim- 
ulations combined with coupled-mode theory. These calculations take fully into account all 
imperfections, including metal losses and imperfect waveguide coupling. We observe that 
there is some optimal R where P is maximized, which corresponds to a balance between 
achieving large coupling between emitter and nano-structure, and ensuring that the plas- 
mon/guide coupling exceeds the enhanced losses at small R. We find that optimal single 
photon efficiencies exceeding 95% are achievable in such a system. 

Such an architecture for quantum communication based on plasmonic devices has several 
important features. First, unlike typical methods of cavity QED, the plasmon excitation 
covers a broad bandwidth and requires no special tuning to achieve resonance. The operation 
speeds can also be quite high because of the sub-wavelength mode volumes associated with 
the plasmons. Finally we note t 



for nanowires 



lally 



29, 3Cj, nanotips 



; rat rapid advances in recent years in fabrication techniques 
I, and sub- wavelength dielectric waveguides J32I, [jjj puts 
such a system in experimental reach. 
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FIG. 1: (a) An emitter coupled to a nanowire is optically excited and decays with high probability 
into the plasmon modes of the nanowire. A single photon source is created by evanescently coupling 
the nanowire to a nearby dielectric waveguide over a length L ex . (b) A similar setup, where a dipole 
emitter is coupled to a metallic nanotip that expands to some final radius R and is then coupled 
to a dielectric waveguide. 




FIG. 2: (a) kn for plasmon modes on a silver nanowire as a function of wire radius R, in units of k\ = 
ujQ^fe\/c. The fundamental plasmon mode (in black) is characterized by a kn<x\/R dependence. 
Inset: the ratio Re fcii/Im kn for the fundamental mode, (b) Solid line: Probability of emission into 
non-plasmon channels, 1 — r p i/(r' + r p i), for a nanowire as a function of R. Dashed line: probability 
of emission into non-plasmon channels vs. v for a nanotip. Dotted line: 1 — r p i(i?)/(r' + r p i) at 
final radius R for a nanotip. Solid points: numerically optimized values of 1 — r p i(i?)/(r' + r p i) for 
a nanotip, obtained via boundary element method simulations. Inset: same plot, zoomed in near 
R,v = 0. 
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FIG. 3: (a) Normalized energy flux for an emitter positioned (from top to bottom) at distances 
hod = 0.002,0.2,0.7 (denoted by the blue circles). The nanotip (whose surface is indicated by 
the dotted lines) has final radius koR = 0.3 and curvature parameter kov = 0.022. The first 
plot is mostly dark and indicates that the emitter decays primarily non-radiatively. The middle 
plot demonstrates efficient excitation of guided plasmons at the final radius R, while the last plot 
exhibits the typical lobe pattern associated with radiative decay, (b) Optimized efficiency of single 
photon generation vs. R. We have assumed that coupling to waveguide modes other than the 
fundamental mode is negligible, i.e., the waveguide is effectively in the single-mode regime. Solid 
line: theoretical efficiency using a nanowire. Dotted line: theoretical efficiency using a nanotip. 
Solid points: efficiency based on numerical simulations of emission near a nanotip, combined with 
coup led- mode equations. 
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